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Abstract 

We present parameter-free tunneling calculations for various cluster radioactivities including 
the diproton decays of atomic nuclei. An uniform folded cluster potential has been suggested 
that is based on a self-consistent mean-field model, with the folding factor determined using the 
quantization conditions of the quasibound cluster state. We have investigated the a-particle and 
heavier-cluster decays of trans- 100 Sn and trans- 208 Pb nuclei, and the observed diproton emission 
from the proton drip-line nucleus 16 Ne, showing the overall reasonable descriptions of cluster ra- 
dioactivities with calculated half-lives agreeing well with experimental data. We have also predicted 
the properties of yet unobserved cluster decays of the exotic nuclei 112 ' 114 Ba, 104 Te and 38 Ti. 

PACS numbers: 23.60.+e, 23.70.+j, 21.10.Tg, 27.60.+j, 27.90. +b 
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The charged-particle emission is one of the most important decay modes of atomic nuclei. 
Glancing at the chart of nuclides, one can find that almost all observed proton-rich exotic 
nuclei starting from A ~ 150 have a radioactivities, a decays occur more widely in nuclei 
heavier than lead isotopes. By a simple picture, the a decay can be understood as a process 
of the a particle being formed in the mother nucleus and emitted tunneling through a 
Coulomb barrier which is created due to the Coulomb interaction between the cluster and 
the remaining nucleons (i.e. the daughter system). The study of a decays has provided rich 
information about the structures of nuclei. 

The emissions of heavier clusters, such as 14 C, 20 O, 24 Ne, 28 Mg and 32 Si, have been well 



established experimenta 



magic nucleus of 208 Pb 



ly in trans-lead nuclei decaying into daughters around the doubly 
. 3, A second island of heavy-cluster radioactivities was pre- 
dicted 0| in trans-tin nuclei decaying into daughters close to the doubly magic nucleus of 
100 Sn. Intense studies, particularly on the most promising case of the 12 C emission from 
the proton drip-line nucleus 114 Ba decaying into 102 Sn, have been made both experimen- 



tally 
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fly 
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12j |. Experiments have not pinned down the 
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re predictions of the partial half-life 
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However, the recent experi- 



and theoretically 
observation of the 12 C decay of 114 Ba. Theoretical 
can be different by several orders of magnitude jj], 
ment |8[ has derived the Q value for the possible 12 C decay of 114 Ba, which is very important 
for the theoretical prediction because the calculated half-life depends dramatically on the Q 
value. At the proton drip line, another exciting phenomenon is the diproton radioactivity of 
nuclei. Modern facilities have been opening new perspectives for the study of exotic decays. 

The decay process of a charged cluster can be treated with the Wentzel-Kramers-Brillouin 
(WKB) approach. In the present work, we have investigated cluster decays from ground 
states to ground states in even-even nuclei. In such decay, the cluster does not carry any 
angular momentum, i.e. L = 0. The decay width can be written as (see, e.g., fis| ) 



h 2 

T = PF— exp 
4/i 



drk(r) 



ri 
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where P is the preformation probability of the cluster being formed in the mother, and fi is 
the reduced mass of the cluster-daughter system. The normalization factor F is given by 
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and the wave number k(r) is defined by 



Hr) = ^\Q-V(r)\, (3) 

where Q is the decay energy and V(r) is the cluster potential in which the cluster moves. r\ 
and r 2 are the classical inner and outer turning points, respectively, obtained by V(r) = Q. 
Then the half-life of the cluster decay is obtained by T1/2 = hln2/T. 

In the above method, the determination of the cluster potential V(r) is a key step for 
the decay calculation. Phenomenologic potential s by Buck et al. [3] and potentials based 



lotentiais Dy 

0HQ 



on effective nucleon-nucleon interactions [14, 15, 16j have been successfully applied to the 
WKB calculations of the charged-particle decays of nuclei. The cluster potential can be 
written as 

V(r) = V N (r)+V c (r), (4) 

where V^(r) and Vc(r) are the nuclear potential and the Coulomb potential, respectively, 
between the cluster and the daughter system. In the present work, the nuclear potential is 
constructed by nucleon potentials with multiplying a folding factor A as follows 

V N (r) = X(N c v n (r) + Z c v p (r)), (5) 

where N c and Z c are the neutron and proton numbers of the cluster, respectively; v n (r) and 
v p (r) are the neutron and proton potentials (excluding the Coulomb potential) respectively, 
obtained from a self-consistent mean-field model. Folding processes have been employed in 



ot 



lovea m 

0,0, 



ier forms of a-nucleus potentials based on effective nucleon-nucleon interactions 

The Coulomb potential Vc(r) is well defined physically and should not be folded. We 
have approximated the Coulomb potential by Vq(t) = Z c v c (r), where v c (r) is the proton 
Coulomb potential obtained from the mean-field calculation. 

In the tunneling model, the Bohr- Sommerf eld quantization condition defines a relation 
between the potential V(r) and decay energy Q for a quasibound cluster state. We have used 
the condition to determine the folding factor A at a given Q value. It will be seen that thus 
determined folded a-nucleus potential leads to a realistic volume integral (per nucleon pair) 
which is obtained from a-particle scatterings experimentally 
condition is given by [3| 



14 , 15j . The Bohr-Sommerfeld 



£ drJpQ - V(r)\ = (2n + 1)| = (G + 1)|, (6) 



where n is the node number in the L = radial wave function of the relative motion of 
the cluster in the potential, and G is the oscillator quantum number with G = 2n 13]. 
The quantum number G can be determined by the Wildermuth condition, G = Ylt=i 9i 15], 
where A c is the mass number of the cluster and gi is the oscillator quantum number of a 
cluster nucleon orbiting in the mean field. Actually, g$ is just the main quantum number in 
the Nilsson labelling. For the lowest-energy cluster decay, cluster nucleons should occupy 
orbits immediately above the Fermi levels of the daughter nucleus. For trans- 100 Sn nuclei, 
for example, cluster nucleons occupy the g = 4 orbits, making G = 4A C . In trans- 208 Pb 
nuclei, the orbits occupied by cluster nucleons have g = 5 for protons and g = 6 for neutrons, 
leading to G = 5Z C +6N C . Hence, the combined use of the Bohr-Sommerfeld and Wildermuth 
quantization conditions can determine the folding factor A. 

Now it has been seen that the present method does not bring any extra adjustable pa- 
rameter. The parameter-free calculation is particularly useful for exotic nuclei where exper- 
imental data are lacking. Calculations based on self-consistent mean-field models should in 
principle give more consistent descriptions on both the structures and decays of nuclei. It has 
been known that mean-field models can in general produce the structure properties of exotic 
nuclei jl^. Mean- field potentials, due to their self-consistent interdependence of proton and 
neutron densities that are fed back into potentials, automatically contain a dependence on 
nucleon numbers, i.e. an isospin dependence, in a self-consistent manner [rs| . 

In a very recent work Q], a form of the a-nucleus potential by V(r) = 2t> p (r) + 2v n (r) 
has been used with v p (r) and v n (r) calculated by mean-field models. However, an extra 
parameter determined by fitting decay data is needed for a-decay calculations • Also, it 
will be seen that the a-nucleus potential without folding gives a too deep well, leading to 
a too large volume integral. Hence, the folding process is also to reduce the depth of the 
cluster potential. 

In our calculations, a preformation factor of P = 1 has been assumed for the various 
cluster decays of even-even nuclei, as suggested in Refs. jll|, [l^. The potentials v n (r), 
v p (r) and v c (r) have been calculated in the spherical daughter system using the Skyrme- 
Hartree-Fock (SHF) approach with the SkI4 force that has been developed with good isospin 
properties Q|. Actually, we found that different Skyrme forces in general lead to similar 
results for cluster decay calculations. Spherical shapes allow us to perform simple one- 
dimension tunneling calculations, as in Refs. HQ- 

Compared with the 
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preformation factor, the Q value is a much more crucial quantity that affects the calculated 
decay half-life, since the half-life is exponentially dependent on the Q value. Therefore, 
we have adopted experimental Q values that can be obtained from the measured masses 
(binding energies) of mothers, daughters and clusters. 

As a test, we have calculated the a-decay property of the spherical nucleus 212 Po that has 
been shown to have a significant a-cluster structure outside the core of the doubly magic 



nucleus 208 Pb 



2l| . Our WKB calculation leads to Tw 2 = 89 ns against the experimental 



half-life of 299 ± 2 ns 



221 ]. The ratio between the calculated and measured half-lives gives 



a preformation factor of 0.3 agreeing with the value of 0.3 calculated by the shell-cluster 
model (2]]]. The folding factor A is determined to be 0.595. Such a folded a-nucleus potential 
gives a rather reasonable volume integral of J v = 325 MeV fm 3 which is in a good agreement 
with the experimental J v ~ 300 — 350 MeV fm 3 for a wide range of nuclei Q], obtained 
from a-particle scattering measurements. The a-nucleus potential without folding [l^ leads 
to a too large volume integral of ~ 550 MeV fm 3 . 

With the development of experimental techniques, more and more exotic nuclei far from 
the stability are being produced. Their decay properties are attracting great interest ex- 
perimentally and theoretically. In the present work, we have investigated the cluster decays 
of trans- 100 Sn and trans- 208 Pb nuclei in which a-particle and heavier-cluster radioactivities 
have been observed or expected. Table |U shows calculated a-decay properties. The calcu- 
lated half-lives agree with experimental data within a factor of ~ 2, except 114 Ba for which 
the experiment of Ref. f| gave a partial half-life of T^ 2 > 1.2 x 10 2 s. As predictions, we 
have also calculated the a-decay half-lives of the unknown proton drip-line nuclei 104 Te and 
112 Ba, shown in Table |U In the calculations, theoretical Q values given by Moller et al. j^J 
have been adopted. However, a change of 1 MeV in the Q value can lead to a variation of 
~ 3 orders of magnitude in the half-life calculation. From Table HJ it can be seen that the 
folding factors of the a-nucleus potentials are quite stable at A ~ 0.5 for trans- 100 Sn nuclei 
and A ~ 0.6 for trans- 208 Pb nuclei. Resulting volume integrals are J v ~ 290 and ~ 330 MeV 
fm 3 for the trans- 100 Sn and trans- 208 Pb nuclei, respectively. 
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Heavier-cluster radioactivities in trans- 208 Pb have been observed well LU, |2|, |3j . Theoret- 
ical investigations have also been made 



24j . Our parameter-free calculations for the 



partial half-lives of the heavy-cluster decays agree with observations within one order of 
magnitude, shown in Table |Hj The determined folding factors have similar values to that 
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TABLE I: Calculated half-lives of the a decays in trans- 100 Sn and trans- 208 Pb nuclei, compared with 
experimental half- lives obtained from evaluated data in Ref. |22|. For 114 Ba, another experimental 

D 



half-life [6j is given for comparison. 
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"Theoretical value from 



of corresponding a-nucleus potentials. For the 12 C radioactivity in 114 Ba, the recent exper- 
iment has derived the decay energy of Q = 19.00 ± 0.04 MeV [8j. With the experimental Q 
value, we have calculated a partial half-life of T\i<i = 1.5 x 10 10 s and a branching ratio of 
b = 4.7 x 10 -8 (b = r( 12 C)/r(a)). The experiment by Guglielmetti et al. [?| has estimated 
a partial half-life of Ti/ 2 > 1.2 x 10 4 s and a branching ratio of b < 3.4 x 10 -5 for the 12 C 
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TABLE II: Partial half-lives and branching ratios (relative to the a decay) of the heavy-cluster 
decays in trans- 208 Pb and trans- 100 Sn nuclei. For trans- 208 Pb nuclei, measured partial half-lives 

n 

from Ref. 



have been compiled in Ref. |llj], and experimental branching ratios and Q values have been obtained 
22]. For trans- 100 Sn nuclei, numbers in square brackets indicate references from which 



the adopted values come. 
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decay of 114 Ba. The nucleus 112 Ba has not been known experimentally. We have used the 
theoretical decay energy of Q = 23.17 MeV given by Moller et al. 23] to calculate the partial 
half-life of the -C decay of -Ba, giving T 1/2 = 67 s. Kumar et al. gave another calculated 
value of Q = 21.46 MeV jlfj|, leading to a half-life of 9.5 x 10 4 s by our calculation. 

Around the proton drip line, the diproton radioactivity is another exciting challenge 



in both experiment and theory. Experiments have observed a few examples of diproton 
emissions from light proton drip-line nuclei 22]. Nazarewicz et al. have investigated 
diproton radioactivities around the doubly magic nucleus 48 Ni in the framework of the WKB 
method with the form of diproton potential by V2 P (r) = 2v p (r) (They included the Coulomb 
potential in v p (r)). The proton potential v p (r) was calculated using various mean- field 
models [20]. They compared the depths of potentials, resulting in a modification on the term 
of \Q — V 2p \ by multiplying an effective mass of m*/m < 1 j20j. Similarly to cluster decays 



discussed above, we have approximated the diproton potential by V2 P (r) = A x2v p (r)+2v c (r) 
(Note that our v p (r) excludes the Coulomb potential). Again, the half-life calculation of the 
diproton decay is dramatically dependent on the decay energy of Q 2p - 

The diproton emission from the proton drip-line nucleus 16 Ne decaying into the proton 
magic nucleus 14 has been observed with an intensity of 100% In this nucleus, the 
two protons occupy the lcfa/2 orbits above the Z = 8 closed shell, leading to G = 4. With 



the experimental decay energy of Q 2p = 1411 ± 20 keV 22], we obtained a folding factor 
of A = 0.662 using the Bohr-Sommerfeld condition. The calculated half-life of the diproton 
emission is 6.5 x 10~ 20 s agreeing with the observed half-life of 9 x 10~ 21 s (3] within one order 
of magnitude. Another promising candidate for the diproton radioactivity is the drip-line 
nucleus 38 Ti decaying into the proton magic nucleus 36 Ca. The experiment has estimated 
a upper limit of 120 ns for the half-life of the diproton decay |25[. With the evaluated 
decay energy of Q 2p = 960 ± 260 keV from Ref . 22] , we have determined a folding factor 
of A = 0.75 at G = 6, resulting in the calculated diproton half-life of T\/ 2 = 120 ns. If the 
evaluated error of ±260 keV in the Q value is considered, our calculations lead to a range 
of T 1/2 = 0.26 ns to 1.7 x 10 6 ns. 

In summary, an uniform folded cluster potential has been suggested for the WKB cal- 
culations of various cluster decays in atomic nuclei. The folded potential is constructed by 
nucleon potentials obtained from the self-consistent Skyrme-Hartree-Fock model, with the 
folding factor determined using the Bohr-Sommerfeld quantization condition combined with 
the Wildermuth condition. This leads to the consistent descriptions of the cluster decay and 
shell structure of a nucleus. The folded a-nucleus potential gives a reasonable volume in- 
tegral, agreeing well with the experimental data obtained from a-scattering measurements. 
No adjustable parameter has been involved in cluster decay calculations, which is particu- 
larly useful for exotic nuclei since their properties are largely unknown experimentally. We 
have investigated the a-particle and heavier-cluster decays of trans- 100 Sn and trans- 208 Pb 
nuclei, and the diproton emission from the proton drip-line nucleus 16 Ne above the magic 
14 0. The calculated half-lives agree with experimental data within a factor of ~ 2 for a 
decays, and within one order of magnitude for heavier-cluster and diproton decays. We have 
also predicted the half-lives of possible a decays in the unknown drip-line nuclei 104 Te and 
112 Ba, the expected 12 C emission from 114 ' 112 Ba, and the diproton decay of 38 Ti. In our 
calculations, deformation effects on cluster radioactivities have not been taken into account, 



8 



which will be discussed in our other works. 
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